function  = Sphere_DB(dX, dY, dZ,...    % 观察点坐标
                      dXc, dYc, dZc,... %球心坐标
                      dR,...       %球体半径
                      dK,...      %磁化率
                      dMr, dIr, dDr,... %剩余磁化强度的大小， 倾角和偏角
                      dA,...    %x轴正方向的方位角
                      dBO, dI, dD)    %背景场的大小， 倾角和偏角 (度)

NT2T = 1E-9;	%
T2NT = 1E9;     %
MU = 4*pi*1E-7; %真空磁导率

% 将角度转换为弧度
dIr = Deg2ArcDeg(dIr);
dDr = Deg2ArcDeg(dDr);
dA = Deg2ArcDeg(dA);
dI = Deg2ArcDeg(dI);
dD = Deg2ArcDeg(dD);

dNx = cos(dI)*cos(dD-dA);
dNy = cos(dIr)*sin(dD-dA);
dNz = sin(dIr);

dLx = cos()*cos()
dLy = 
dLz = 

dBO = dBO*NT2T;
dMi = dK*(dBO/MU);

dMx = dMi*dNx + dMr*dLx;
dMy = dMi*dNy + dMr*dLy;
dMz = dMi*dNz + dMr*dLz

dVxx = Vaa(dX, dY, dZ, dXc, dYc, dZc);
dVyy = Vaa(dY, dZ, dX, dYc, dZc, dXc);
dVzz = Vaa(dZ, dX, dY, dZc, dXc, dYc);
dVxy = Vaa(dX, dY, dZ, dXc, dYc, dZc);
dVyz = Vaa(dY, dZ, dX, dYc, dZc, dXc);
dVzx = Vaa(dZ, dX, dY, dZc, dXc, dYc);

dBOx = dBOx+dBax; 
dBOy = 
dBOz = 

dBx = dBOx+dBax;
dBy = dBOy+dBay;
dBz = dBOz+dBaz;

dDB1 = T2Nt * sqrt(dBx^2 + dBy^2 + dBz^2);
dDB2 = T2NT * dBax*dNx + dBay*dNy+dBaz*dNz;

end

